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ABSTBACT 


The  problems  in  megiietohydrodynamlce ,  in  which  an  inviscid  compressible 
fluid  with  small  conductivity  flows  steadily  past  a  slender  body  of  an  arbitrary 
cross  section,  in  the  presence  of  an  applied  magnetic  field  parallel  to  the  uni** 
form  flow  are  considered.  By  the  use  of  the  slender  body  approximation  we  shall 
discuss  the  character  of  the  velocity  and  the  magnetic  fields,  and  obtain  the  drag 
and  the  lateral  force  exerted  on  the  body.  For  examples,  the  flows  past  a  body  of 
an  elliptic  cross  section  and  a  body  of  revolution  at  a  small  incidence  will  be 
discussed.  It  is  noticeable  that  the  induced  drag  for  a  slender  body  of  revolu¬ 
tion  with  pointed  nose  and  tall  ends  is  negative,  and  the  lift  for  the  same  body 
takes  the  negative  sign  for  a  positive  incidence. 
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INTRODUCTION 


Recently  Sears  and  Resler  (1)  investigated  the  flow  of  a  non-vlscous  in¬ 
compressible  fluid  with  an  electrical  conductivity  past  a  two-dimensional  thin  ob¬ 
stacle  by  the  linearized  theory.  In  their  paper,  the  conductivity  of  fluid  is 
supposed  to  be  infinite  or  very  high.  The  same  problem  for  an  arbitrary  conduc¬ 
tivity  was  studied  by  McCune  (2).  The  correlation  effect  between  the  compressi¬ 
bility  and  the  conductivity  was  studied  by  Taniuchl  (3),  Kogan  (4),  McCune  and 
Resler  (5),  and  Sears  (6).  In  these  papers,  the  conductivity  of  fluid  is  supposed 
to  be  infinite.  The  flows  of  the  non-viscous  compressible  fluid  with  a  small  con¬ 
ductivity  past  a  two-dimensional  thin  body,  and  a  body  of  revolution  at  zero  inci¬ 
dence  were  studied  by  Sakural  (7),  (8),  and  Kusukawa  (9),  respectively.  Sakural 
and  Kusukawa  take  the  magnetohydrodynamlc  Stokes'  and  Oseen's  approximations,  re¬ 
spectively.  Resler  and  McCune  (10)  studied  the  flow  of  a  compressible  fluid  with 
an  arbitrary  conductivity  past  a  sinusoidal  wall.  Ando  (11),  (12),  discussed  the 
character  of  the  linearized  equation  for  the  flow  of  a  compressible  fluid  with  an 
arbitrary  conductivity  past  a  thin  body,  in  the  presence  of  the  uniform  applied 
magnetic  field  in  an  arbitrary  direction. 

In  the  present  paper  we  shall  study  the  flow  of  a  non-viscous  compressible 
fluid  with  a  low  conductivity  past  a  slender  body  of  an  arbitrary  cross  section, 
in  the  presence  of  a  uniform  applied  magnetic  field  parallel  to  the  uniform  flow. 

We  shall  linearize  the  fundamental  equation,  and  fpply  the  slender  body  theory 
(13)  to  the  present  magnetohydrodynamlc  flow.  The  induced  electric  field  is  quasi- 
two-dimensional  in  the  sectional  plane.  Both  the  velocity  and  the  magnetic  fields 
can  be  decomposed  into  the  axial  field  and  the  irrotational  two-dimensional  field 
in  the  sectional  plane.  He  shall  discuss  the  circulation  of  fluid.  The  total 
force  exerted  on  the  body  will  be  calculated.  For  an  example,  the  drag  and  the 
lift  for  a  body  of  revolution  with  pointed  nose  and  tall  ends  at  a  small  incidence 
will  be  studied.  It  is  worthwhile  to  notice  here  that  the  induced  drag  is  nega¬ 
tive,  and  the  lift  is  negative  for  a  positive  incidence. 
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I  FUNDAMENTAL  EQUATI(»IS  FOR  FLUID  FLOW 


We  shall  consider  the  flow  of  a  non-viscous  compressible  fluid  with  a  low 
electrical  conductivity past  a  body  of  an  arbitrary  cross  section,  In  the 
presence  of  an  applied  uniform  magnetic  field  H*  parallel  to  the  uniform  flow 
\J*  .  We  shall  take  the  cms  electromagnetic  system  of  units  and  the  cartesian 

coordinates  x*,  y*,  s*,  where  s*  Is  parallel  to  the  uniform  flow.  For  simplicity, 
we  shall  normalise  the  physical  quantities  with  respect  to  the  length  L*  of  the 
body,  the  uniform  velocity  U*  >  the  uniform  magnetic  field  H*  ,  and  the  density 
and  the  temperature  ~J*  In  the  uniform  flow  such  that 


X« 


1. 


L* 


3 


2J! 


IT-  ^ 

TT 

a= 

u* 

u: 

H* 

“  HJU' 

■ur* 


H?  u: 


C, 


A  =  -fr  ii. 

H* 


L* 


h: 


where  u,  v,  w,  denote  the  components  of  the  perturbation  velocity,  iia  ,  ,K» 

the  perturbation  magnetic  field, ew  ,  Cy  ,  the  Induced  electric  field,  -jn  ,  i't 

,  the  Induced  electric  current  density,  p  the  pressure,  (»  the  density,  T  the 
temperature,  and  the  notations  with  and  without  asterisk  represent  the  physical 
and  the  dimensionless  quantities,  respectively.  For  convenience  we  shall  Intro¬ 
duce  the  dimensionless  parameters  defined  by 

H.** 


S  = 


4ir 


SR-.- 


.•*un* 


where  S  and  Rm  are  the  pressure  number  and  the  magnetic  Reynolds  number,  respec¬ 
tively. 

Supposing  that  the  magnetic  permeabilities  of  both  the  body  and  the  fluid 
are  unity,  and  the  fluid  Is  non-vlscous  and  thermally  non-conducting,  we  can  write 
the  equation  of  motion,  the  equation  of  continuity,  the  Ohm's  law,  the  Maxwell's 
equation,  and  the  conservation  of  energy  In  the  following  form; 
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H = f  1'*-’ 

-  IS  --  if  *  5?-  =  f  h>-  -  j' 

5:<f“->'-|r(f'"^*li(f’H  •° 

i,  •  R7»^v(l*h,)-(ltwH,  »«x] 
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ao 
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(4a) 


(4c) 


(4e) 


(4g) 


ds  ~  ax 
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_ 

^ 

as  ” 


o>g» 
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ac* 
d  s 


a»  O 


(la) 


(lb) 


(Ic) 


(2) 


(3a) 


(3b) 


(3c) 


(4b) 


(4d) 


(4f) 


(4h) 


(5) 


d«x  _  3  C* 

where  Cv  denotes  the  specific  heat  of  tl\e  fluid. 

When  the  conductivity  O*  of  the  fluid  is  low,  the  coupling  between  the 
velocity  and  the  magnetic  field  is  small.  Actually  the  uniform  applied  magnetic 
field  is  scarcely  distorted  for  a  very  small  conductivity.  As  the  conductivity 
increases,  the  coupling  becomes  strong.  It  is  well  known  that  the  magnetic 
field  is  frozen  into  the  material  in  the  extreme  case  of  infinite  conductivity. 
From  these  considerations,  we  can  suppose  that  components  of  perturbation  magnetic 
field  ,  /ly  ,  Ut  are  the  order  of  magnitude  of  the  thickness  ratio  *2^,  of 

the  body,  at  most,  for  a  not  so  large  conductivity , such  that 

for  Q  £  ]  anil  Rwi  »  I 

0(l.x)iO(t.)  0(^)i6(T.)  oCk)iO(Z:')  (6a) 


On  the  other  hand,  we  can  estimate  the  order  of  magnitude  of  the  components 


3 


of  the  perturbation  velocities  and  the  density  near  the  surface  of  the  body  based 
on  the  knowledge  of  the  conventional  flow  past  a  slender  body  without  a  magnetic 
field  (13)  such  that 

o('u)  =  o(2:)  oM^oCx)  o(ut)  =  o('z;*Io3  4) 

,  ‘  (6b) 

O(^)  = 

Remembering  these  order  estimations,  and  retaining  the  leading  terms,  from  Equa¬ 
tions  (3a),  (3b)  and  (3c)  we  obtain 


Rk.  f  kx 

-js  = 


(7a) 

(7b) 

(7c) 


Since  every  term  In  each  equation  may  be  considered  to  be  the  same  order  of  magni¬ 
tude,  we  have 

o(e.)xo(r.')  o(€>,)-:OC%)  o(e,^  =  o(r.O 


We  shall  take  a  line  element  u  along  a  stream  line  defined  by 


(6c) 


u. 


dx. 

dX 


2L  djL 
%  ~ 


i-hlM- 


ds 

dl 


where  f  denotes  the  magnitude  of  the  nondlmenslonal  velocity.  Integrating  Equa¬ 
tions  (la)-(lc)  with  respect  to  €  along  a  stream  line,  and  remembering  the  order 
of  estimations  (6a)-(6c)  we  shall  obtain 

=  -  Zw-(u*- v')  -2  «jr{('»r-Ay)(v'-A,  +  Cx)+(k»*ta,)(Kx-u+ey)jd^  (8a) 

where  denotes  the  pressure  coefficient  defined  by 

_  -r*-  T.* 


since  the  body  is  very  thin,  and  stream  lines  are  almost  parallel  to  the  uni¬ 
form  flow  velocity,  we  can  approximately  replace  the  variable  .Z  by  S.  This  is  the 
quadratic  formula  for  the  pressure  coefficient.  If  we  put  A.«o  in  Equation  (8a) 
we  have  the  pressure  formula  for  the  conventional  case  (14).  By  the  use  of  Equa¬ 
tion  (8a)  we  can  calculate  the  pressure  distribution,  if  we  find  the  velocity  and 
the  magnetic  field  over  the  body. 

Especially  in  the  axially  symmetric  case,  the  Induced  electric  field  can  be 
shown  to  be  zero  because  of  the  geometrical  synsnetry . (9)  Thus  Equation  (8a)  Is 
reduced  to 


cr  =  -2vr-(u%»-'j-  || 


ds 


(8b) 
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The  last  teim  of  the  right  hand  side  of  this  equation  represents  the  energy  loss 
due  to  Joule  heat.  The  pressure  drop  is  connected  with  the  irreversible  process. 


II  LINEARIZATION  OF  EQUATIONS  FOR  FLUID 


We  shall  linearize  the  fundamental  equations  for  the  fluid  flow  obtained  in 
the  last  section.  If  we  neglect  the  terms  involving  the  square  of  perturbation 
quantities,  the  equations  of  motion  and  the  equation  of  continuity  are  reduced  to 


di  Ts 


2if: 

Ps 

d  s 


(9a) 

(9b) 

(9c) 

(10) 

(lla) 

(llb) 


Combining  Equation  (4a)  with  Equation  (7a),  we  obtain 

Similarly,  from  Equations  (4b)  and  (7b),  we  obtain 

Considering  the  order  estimations  (6c),  we  shall  find  that  the  axial  components  of 
the  electric  fields*  and  the  electric  current  density can  be  safely  neglected 
in  comparison  with  the  x-  and  the  y  components  within  the  accuracy  of  the  linear¬ 
ized  theory.  Thus  (4c)  and  (4h)  can  be  respectively,  approximated  by 

Equations  (4f)  and  (4g)  imply  that  and^*/^»  vanish  in  the  present  approxi¬ 

mation.  Thus  the  Induced  electric  field  forms  a  quasi-two-dimensional  field.  Re¬ 
membering  Equations  (6c),  we  shall  find  that  the  last  parenthesis  of  the  right 
hand  side  of  Equation  (5)  is  the  order  of  o(flLir,*  ).  Neglecting  the  higher  order 
terms,  we  have  m.-r*  y 

***  p  =  -X  — — 

Combining  this  equation  with  the  equation  of  state  of  an  ideal  gas,  we  have 


fV 


const 


where T^denotes  the  ratio  of  the  specific  heats  Vcv  .  The  fluid  conforms  the 
law  of  the  adiabatic  change  within  the  accuracy  of  the  linearized  theory. 
Introducing  the  Mach  number  of  the  uniform  flow 

^ 

where  is  the  free  stream  speed  of  sound,  we  have 

(12) 

From  Equation  (11c)  we  shall  introduce  the  quasi-two-dimenslonal  magnetic 
potential  such  as 
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*»y=  3'f  (13a) 

Eliminating  from  Equations  (11a)  and  (11b),  and  considering  Equations  (11c)  and 

3v  3-u 


(lid),  we  obtain 


ax 


■=  o 


From  this  equation,  we  can  Introduce  the  quasi-two-dlmenslonal  velocity  potential 
defined  by 

3  IP  _  3  If 

ax  5^  (13b) 


From  Equation  (lid),  we  shall  Introduce  the  electric  force  function  such  as 

3».0<  _  3  vfl« 

=  =  (13c) 

Combining  these  equations  described  with  each  other,  we  shall  obtain  the  fundamen¬ 
tal  equations  governing  (f  .  and  i/X  such  that 


a[  { a-  B*  -  B*  ft  It]  fe}  B*  a  J  ^ = 


A  JX  =  o 

and 


(lA) 

(15) 


«^th  A-  ^ 

If  the  uniform  flow  Is  subsonic,  takes  the  negative  sign.  Besides,  the 
unknown  functions  'f  ,  and  Jh  have  to  satisfy  the  following  supplementary  equa¬ 
tions. 


a  s 


The  axial  components  of  the  perturbation  velocity  and  the  magnetic  fields 
can  be  given  by 

(17a) 


=  o 

(16a) 

-o 

(16b) 

<  a  O 

(16c) 

o 

(16d) 

(17b) 


Considering  Equations  (13b)  and  (17a)  we  find  that  the  velocity  field  Is  Irrota- 
tlonal  In  the  sectional  plane  perpendicular  to  the  s-axis,  but  it  Is  rotational  In 
the  meridian  plane.  The  situation  Is  the  same  for  the  magnetic  field. 


Ill  FUNDAMENTAL  EQUATIONS  FOR  INSIDE  ELECTROMAGNETIC  FIELD 
In  the  Interior  of  the  body,  the  Ohm's  law  Is  simply  represented  by 
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=  (18a)  dsb) 

j,yi  =Rie3i 
with  Rw=4‘ir(rt,*U.''Lt 


where^fc  denotes  the  electric  conductivity  of  the  body,  and  the  subscript  "i"  de¬ 
notes  the  quantities  in  the  Interior  of  the  body.  The  Maxwell's  equations  are  rep¬ 
resented  by  Equations  (4a)-(4h).  As  described  above  (c.f.  equation  (6c)),  the 
axial  component  of  the  outside  electric  field  6s  is  negligibly  small.  If  we  con¬ 
sider  that  the  s-component  of  the  electric  field  should  be  continuous  on  the  sur¬ 
face  of  the  body,  and  the  thickness  of  the  body  is  the  order  of  ,  we  can 

conclude  the  s-component 6s«  of  the  inside  electric  field  to  be  the  order  of  <>(%*). 
Remembering  Equation  (18c)  we  can  safely  neglect  6s;  and  within  the  accuracy  of 
the  linearized  theory.  Considering  Equations  (4c)  and  (4h) ,  and  neglecting  the 
higher  order  terms,  we  can  introduce  the  inside  magnetic  potential  and  the  in¬ 
side  electric  force  function  defined  by 


A 

'  alT' 

V-5Y- 

***  ^  Y 

treatment  as 

that  for  the  outer 

— - o  -  - - - - - - - j  ’ —  -  — 

lowing  fundamental  and  supplementary  equations  governing  4^^andj^., 


(19a) 
(19b) 
i  the  i 

(20) 


The  axial  component  of  the  magnetic  field  can  be  given  bj 


(21) 

(22) 


(23) 


IV  GENERAL  SOLUTIONS  OF  FUNDAMENTAL  EQUATIONS 
(A)  OUTER  SOLUTIONS 


First  we  shall  consider  the  supersonic  case.  Since  the  fundamental  equation 
(14)  is  linear,  the  general  solution  should  be  represented  by  the  sum  of  three  in¬ 
dependent  solutions  and  such  that 


The  solution 


satisfies  the  Laplace's  equation. 


(24) 


(25) 

Introducing  the  Laplace  transform  of  the  solution  ,  and  the  Fourier  trans¬ 

form  of  the  solution  with  respect  to  S  ,  we  shall  obtain  the  equations 
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governing 

and  such  that 

(A-v)>?.  =  o 

(26a) 

(a-f)  <Pt=o 

(26b) 

with 

(27a) 

(27b) 

4 

k  \J  ( BN-i)'  cu’  4  2  i  u»  (B'-k0(8’ft-  Rh.')  -  (B‘a+ 

where  ^  and  (a>  denote  the  parameters  of  the  Laplace  and  the  Fourier  transforms, 
respectively.  The  solution  of  Equation  (26a)  satisfying  the  condition  that  the 
perturbation  velocity  should  vanish  at  infinity,  can  be  expressed  such  as 


^ =  2  /»H (X)  K.  {yCr  rhos (»e  4  < ( \)) 


(28) 


where  f^,  ®  ,S  denote  the  cylindrical  coordinates,  and  Vh^A)  the  arbitrary 

functions  of  \  ,  andKt.(V5  r)  Is  the  n-th  order  modified  Bessel  function  of  the 
second  kind.  Near  the  body,  where  r*  is  small,  the  Bessel  functions  can  be  expanded 
into  the  following  series, 

^  (4Vw^r)+ C  - 

-  1-^  + . 

it  . 

^  r^  +  C  + . 

Kh  ( i  [( •^ -•) •  -  (►•  -2) ?  {jTr)  ^  ^“2~  j  n  i3 


with  C  =:  0.5772 

Remembering  (27a),  we  shall  find  that«<^  takes  a  moderate  value  for  a  moderate 
value  of  \  ,  when  the  parameters  B*  ,  Q.  and  Rin  are  the  order  of  unity,  0(1) 

In  such  a  case  we  may  take  the  first  three  terms  as  an  approximation  of  the  Bessel 
functions  for  a  small  value  of  P  .  Moreover,  the  second  and  the  third  terms  in 
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the  series  expansion  of  1C*  are  much  smaller  than  the  first  term.  Since  the  second 
and  the  third  terms  are  not  essential,  they  may  be  safely  replaced  by 


From  similar  considerations,  K,  and  K**  can  be  approximately  represented  by 
Kf,  (TfW  p)  •  Uj  r-  +  (Xm  ‘ 


w*  (Q  O  * 


For  convenience,  we  shall  introduce  the  complex  variable  in  the  sectional  plane 
such  as  •  ■ a 

H  *  3  re" 

Using  the  approximate  expressions  of  I^k  ,  and  applying  the  Inverse  Laplace 
transformation  to  obtain  the  solution  In  the  following  form, 

<t,.  -cP., 


with 


•fai  =  lo^  *■  53  J 

^  a*  r^* 

*•  J*  -<►*■ 


(29a) 


>*9 

i«r  e 


-  or,  **5  *  ^ 


>**3 


a) 


1 


where  and  <**•! (a)  denote  arbitrary  pure  real  and  complex  functions  of  S  ,  R« 

denotes  the  real  part  of  complex  function,  and«^  means  the  linear  operator  corres¬ 
ponding  to  the  expression  (27a).  Similarly,  by  the  use  of  the  Fourier  transforma¬ 
tion,  we  can  obtain  the  solution  In  the  following  form 


(29b) 


The  functions  4**., ‘It,  ,  ,  and^ltg  are  obtained  by  replacing*^  ,  <1*  and  in 

Equation  (29a)  by  ^  , b<  and  fc,.  ,  where ^  means  the  linear  operator  correspond¬ 

ing  to  the  expression  (27b),  and  b.  and  lb»  denote  the  arbitrary  pure  real  and  the 
complex  functions  of  s  different  froma./^)  and  The  general  solution ‘fj  of 

Equation  (25)  satisfying  the  condition  at  Infinity  should  be  given  by 

<^3= [-^  r  ‘‘’3^ 

where  and denote  the  arbitrary  pure  real  and  the  complex  functions  of  s, 
respectively.  Substituting  Equations  (29a),  (29b)  and  (30)  Into  Equation  (24),  we 
shall  obtain  the  solution  such  that 


(30) 
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with 

V.-  'f..*  ‘fi. 

Considering  Equation  (29a),  we  can  find  the  following  relations, 

=  A  ■*  s  O 

and  *  OOty  k,)oC^*) 

oM  =  0lr,^)0(^) 

oC*fg)  ‘  oC't,')  oC^f.) 

Since  the  fundamental  equation  governing  ^  is  the  same  as  that  governing  df 
we  can  obtain  the  solution  ^  in  the  similar  form, 

V.  -  ^>C.+ (32) 

+  +M'Js 

The  functions  +1.  , +4.  ,V'd*  , '4»e.  ,'4»di  ,  ,  ^ci  ,  and  can  be  ob¬ 
tained  by  replacing  the  arbitrary  functions  P,  ,  a,  ,  Olh  ,  ,  and  lb»  in 

the  solution  with  the  different  arbitrary  functions  Pa  ,  Co  ,  Ct,  ,  <1,  ,  and 
do  ,  respectively.  We  have  the  similar  relations  to  (31b)  and  (31c). 

The  solution  of  Equation  (15)  satisfying  the  condition  at  infinity  can  be 
readily  written  in  the  following  form, 

JX  *  “R*  [-  i-  fi  loi  «  +  Z  IKoS  *'1  (33) 

I'd  -* 

Next,  in  the  subsonic  case,  both  the  independent  solutions  and  i  are 
obtained  by  the  Fourier  transformation.  The  results  are  given  in  the  same  expres¬ 
sion  as  those  of  the  supersonic  case  by  replacing  the  linear  operators  ^  and  A 
withY  and  / 

(B)  INNER  SOLUTIONS 

The  solution  of  Equation  (20)  without  any  singularity  along  the  axis 
except  the  logarithmic  singularity  is  given  by 

(34) 

The  solution  M’n*  is  the  harmonic  function  such  as 

VVk,  -  “Ra  f-  t  loe2  S  (35) 

The  Fourier  transform of  the  function*f^  with  respect  to  S  should  satisfy  the 

equation  .  _■«  — 

(n-  U)  V  U»£  s  o 

'  (36) 

The  general  solution  of  Equation  (36)  without  any  singularity  along  the  axis  is  ob 
tained  in  the  following  form 


^•P. -o 

^  A  tPjj  =  o 


(31b) 


(31c) 


(31«) 
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£*  FhM  cos(h®)a»7^(u»)) 

»**  O  V  * 


(37) 


where and^tt(w)  are  the  arbitrary  functions  of  ,  and  ^•.(^^'Odenotes  the  n-th 
order  Bessel  function  of  the  first  kind.  Since  the  thickness  of  the  body  is  very 
small,  the  Bessel  functions  can  be  approximated  by  the  first  two-terms  of  their  ex¬ 
pansions  into  the  r  -power  series,  such  that 

X  (tCP'') 


n !  2"  (*»♦!)•  2'’** 

Using  this  approximation,  and  applying  the  inverse  Fourier  transformation 
we  shall  obtain  the  solution  V)r;  such  that 


tovf, 


with 


'Vi.-Tur 


cs) 


^  1^1 

Considering  the  solution  (38a),  we  shall  find  the  following  relations. 


as’ 


and 


oCH'O  =  o(r.')o('Vf.) 


(38a) 


(38b) 

(38c) 


The  solution of  Equation  (21)  can  be  easily  found  such  that 

**  H 

J  (39) 

As  we  shall  discuss  later,  the  logarithmic  singularity  of  the  inside  solu¬ 
tions  is  caused  by  the  current  flowing  in  an  infinitesimally  thin  filament  of  in¬ 
finite  conductivity  Inserted  along  the  axis  of  the  body. 


V  SUPPLEMENTARY  EQUATIONS 


The  general  solutions  obtained  above  contain  several  arbitrary  functions  of 
S  .  They  should  be  determined  by  satisfying  the  boundary  conditions  and  the  sup¬ 
plementary  equations  (16a)-(16d)  and  (22). 

Remembering  Equations  (31a)-(31c),  and  neglecting  the  higher  order  terms 
•(t;;  we  can  write  Equation  (16a)  such  that 

Substituting  the  expressions  (29a)  into  Equation  (40)  and  comparing  the  co¬ 
efficients  of  the  terms  involving  the  same  function  of  2  >  ve  shall  obtain 


\ 


cl  S 


(41a) 

(41b) 
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and 


Vat  -  8'  [«:  1«) 

Integrating  Equation  (Ala),  and  considering  the  condition  that  all  the  perturba¬ 
tions  should  vanish  at  infinity  upstream,  we  have 

^■^a(r.-ri  +  r3)  =  o 

Similarly,  from  Equation  (16b)  we  obtain 

^tRH.(r.-r^^r3)*o 

and  cIh  ^  D*^J'(s)-h  d/J(s)4-lh:rs)  +i?H.  [aC  (5)  bl  (s) 

By  the  use  of  the  same  procedure.  Equations  (16c)  and  (16d)  are  respectively,  re¬ 
duced  to 

-^<Ih  =  o 

*  pel'  (s)+^X(s)4l?K.  -^dllJJ^O 

)b«^i.  c.  sCo  dl.sd. 

The  supplementary  condition  (22)  for  the  inside  field  supplies 

^  r?  j. 

VI  BOUNDARY  CONDITIONS 

We  shall  take  a  plane  perpendicular  to  the  ^  -axis.  The  cross  section  of  the 
body  in  the  plane  forms  a  closed  curve,  C,  say.  The  outward  directing  normal  vec¬ 
tor  and  the  counterclockwise  tangential  vector  of  the  curve  C  are  denoted  by  in 
and t'  ,  respectively.  Within  the  accuracy  of  the  linearized  theory,  the  boundary 
conditions  can  be  given  in  the  following  form 

(i)  Since  the  fluid  must  flow  past  the  body,  we  have 

(42a) 

If  the  body  is  given,  the  right  hand  side  of  this  equation  is  a  known  function, 
integrating  Equation  (A2a)  along  the  curve  C  with  respect  to  T,  we  obtain 

a.  ( y?  4  b,  Cs)  =  ^  Z '( (^2‘>) 

where  2(s)'^enotes  the  area  of  tie  cross  section  of  the  body. 

(il)  Since  the  magnetic  field  should  be  continuous  on  the  surface  of  the 


(Alg) 

(Alh) 


(Ali) 

(Alj) 


(Ala') 

(Aid) 

(Ale) 

(Alf) 


(Ale) 
4-r  n  >  I 
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(43a) 


surface  of  the  body,  ve  have 

3V  QZ  QZ 

A3.  .  a^t  ^A3±1 

^  d  n  ^  r> 

-R»,(>^’,-'^>.-J^)=4±b.  »  3;tii  +  Rw Ji,, 

ds  as 

Integrating  Equations  (43a)  and  (43b)  along  the  curve  C  with  respect  to  'IS  ,  we 
obtain,  respectively 

Q  =  Hf 


(43b) 

(43c) 


and 


3ir 


G(s) 


(43d) 

(43e) 


G(^  =  j 


When  3'r/jT'  is  much  smaller  than  unity, S(*)  can  be  approximated  by  the  sectional 
area^fs)  of  the  body. 


(lii)  Since  the  tangential  component  of  the  electric  field  should  be  contln 
uous  on  the  surface  of  the  body,  we  have 

DJt  ^ 

Q  h  Ah 


(44) 


induced  on  the  surface  of  the  body 


The  total  charge  surface  density  (SJ* 
should  be  given  by  g-  _  AJi-i 

'  "  (45a) 

with  nr  - 

’  H,*  U 

Integrating  Equation  (45a)  along  the  curve  C  with  respect  to ,  we  shall  obtain 

^  cl^  =•  27r  (fj- r^)  (45b) 

The  left-hand  side  of  this  equation  denotes  the  amount  of  the  surface  charge  per 
unit  width  in  the  s-directlon. 

The  total  amount  of  the  surface  charge  can  be  evaluated  by 


with 


o;  d  r  d  s  =  2ir  Jjfj  (5)  -  (5)  J  a  s 

4ir  ^ 


(45c) 


1'  = 


HTU/U* 

(iv)  We  shall  consider  the  condition  at  the  axis  (r=o). 

(a)  In  the  usual  case  we  have  no  singularity  within  the  body.  All  the  vari 
ables  should  be  continuous  and  finite  in  the  interior  of  the  body. 

(b)  We  shall  consider  such  a  geometrical  configuration  that  an  infinitesi¬ 

mally  thin  filament  with  infinite  conductivity  is  inserted  along  the  s-axls.  The 
electric  current  flows  in  the  filament  without  any  electric  field  along  the  fila¬ 
ment.  Moreover  some  electric  charge  can  be  induced  along  the  filament,  and  the 
source  singularity  of  the  electric  field  at  can  be  permitted.  As  a  result, 

we  may  have  a  circular  magnetic  field  and  a  circulation  of  the  fluid. 

Considering  the  continuity  of  current,  we  have 

j=2irr, 
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*T-* 

where  denotes  the  electric  current  flowing  in  the  filament.  The  length  density 
Pi  of  the  electric  charge  induced  along  the  filament  can  be  obtained  by 

P,- =  2-rrfV  (47) 

Remembering  the  fact  that  f*  and  Q  denote  the  source  strengths  of  the  inner  and 
the  outer  radial  electric  fields,  respectively,  we  car  realize  the  physical  mean¬ 
ing  of  Equations  (45b)  and  (47).  For  example,  we  shall  consider  the  case,  where 
both  ri  and  are  positive.  In  the  interior  of  the  body  lines  of  electric  force 
corresponding  to  the  strength  Ft  starts  from  the  filament  and  terminates  at  the 
surface  of  the  body.  In  the  outside  of  the  body  lines  of  electric  force  corres¬ 
ponding  to  the  strength  starts  from  the  surface  toward  infinity. 


VII  CIRCULATION 

We  shall  obtain  the  five  relations  among  the  functions  ^  ,  Pi  ,  C  .U 
and  I*-  represented  by  Equations  (41a'),  (41d),  (411),  (43d)  and  (46). 

(a)  If  there  is  no  filament  along  the  axis,  we  can  put 

J(5)  (^8) 

Combining  Equation  (48)  with  the  five  relations  described  above,  we  shall  obtain 

(49) 

Thus  we  have  found  that  the  circulation  of  the  fluid  vanishes  together  with  the 
singularity  at  the  axis. 

(b)  We  shall  consider  the  case,  where  the  filament  is  inserted  along  the 
s-axis,  and  terminates  at  the  both  pointed  ends  of  the  body.  In  this  case  we  can 
put 


J(5i)  =0  for  S  £  O 
J(S)  #  O  for  O  <  s  <  / 


and  I  s  S' 

(50) 


where  physically  considering, should  be  a  continuous  function  of  S  .  Combining 
Equation  (50)  with  the  five  relations  amoung  P  .  >  and  ,  we 

obtain 


(51a) 

i;=T  ^ 

(51c) 

j-f  r 

(51e) 

in 

d  s 


a=»SRm 


ds 


(51b) 
(5  Id) 


By  the  use  of  Equations  (45c)  and  (47)  we  can  calculate  the  total  amount  of 
the  electric  charge  on  the  surface  of  the  body,  and  the  length  density  of  the 
charge  along  the  filament,  respectively,  such  as 

‘j.,  =2Tr  [% -or;-5  ^ as  (5if) 


and 


zw  d n 

ds 


(51g) 
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Thus,  if  the  circulation  T*  of  the  flow  is  determined,  all  other  quantities!^  , [5 » 
[4  ,  rj  ,  J  ,  ,  and  P  can  be  obtained.  It  is  of  interest  to  consider  the 

following  extreme  cases. 

(i)  If  the  conductivity  Cfc  of  the  body  vanishes,  all  r*i  except  should 
vanish.  The  velocity  field  has  no  circulation  as  well  as  the  magnetic  field.  No 
current  flows  in  the  filament.  Remembering  Equation  (39),  we  realize  that  1^5- 
represents  the  radial  electric  field  in  the  interior  of  the  body.  Then  the  inside 
radial  electric  field  may  remain  in  this  case.  As  a  result  and  may  exist. 

(ii)  In  the  case,  where  the  conductivity  S' ^  of  the  fluid  vanishes,  we  find 
that  all  Ps  except  P 3  vanish.  The  velocity  and  the  magnetic  fields  have  no  cir¬ 
culation.  This  case  corresponds  to  the  conventional  case.  The  outside  radial 
electric  field  corresponding  to  Pj  may  remain. 


(iii)  In  a  hypersonic  flow,  in  the  presence  of  an  applied  magnetic  field,  we 
use  the  magnetohydrodynamic  Stokes  approximation  frequently.  In  this  approximation 
we  suppose  that  R^,  andRw  vanish,  but  Q  and  Q 4  still  remain  finite,  because  $  be¬ 
comes  infinite.  The  circulations  and  Hi.  of  the  outer  and  the  inner  magnetic 
fields  vanish,  but  the  circulation  Jl  of  the  flow  and  the  electric  radial  fields 
Ij  and  Pj-  may  remain.  Corresponding  to  this  fact,  the  current  J  vanishes,  while 
the  induced  charge  distributions  and  J?,  remain. 


We  shall  return  to  the  general  case,  where  and  Q.  are  the  order  of  unity 
at  most.  Remembering  Equations  (31b)  and  (42a),  we  realize  that  should  satisfy 
the  two-dimensional  Laplace's  equation  and  the  Neuman’s  boundary  condition.  It  is 
well  known  that  the  solution  Lf,  can  be  determined  except  an  arbitrary  function 

r»  If  the  curve  C  is  smooth,  we  have  no  mathematical  reason  to  determine 
the  value  of  ^(4) .  If  the  curve  C  has  a  sharp  edge , ® /J T become  infinite  at 
the  edge.  We  can  choose  an  appropriate  value  of  P  in  order  to  prevent  the  infi¬ 
nite  velocity.  These  circumstances  are  similar  to  the  conventional  stationary  flow 
past  a  two-dimensional  obstacle  with  a  sharp  trailing  edge.  If  the  fluid  could 
flow  around  the  sharp  corner  in  the  initial  stage  of  the  magnetohydramic  flow,  the 
tangential  velocity  would  take  a  very  large  value.  Then  the  normal  component  of 
the  induced  current  might  become  tremendously  large,  unless  the  body  is  an  insul¬ 
ator,  According  to  the  interaction  between  the  induced  current  and  the  applied 
uniform  magnetic  field,  a  large  Lorentz  force  exerts  on  the  fluid.  Since  the 
fluid  is  detained  by  the  Lorentz  force,  it  cannot  flow  around  the  sharp  edge.  Thus 
there  occurs  an  appropriate  circulation  around  the  curve  C  which  will  make  the  po¬ 
sition  at  the  edge  coincide  with  a  stagnation  point. 


By  the  use  of  the  Fourier  theorem,  we  can  obtain  the  following  result.  The 
outside  region  of  the  curve  0  in  the  z-plane  can  be  conformally  mapped  into  the 
outside  region  of  the  unit  circle^’®  in  the  J  -plane  by  an  appropriate  analytic 
function  f(»)  .  The  function  ids  in  Equation  (42a)  is  expressed  as  the  func¬ 
tion  of  &  and  S  such  that 


ds 


Ef«-  5) 


(52a) 


If  the  functions  Hi)  and  e(e.')  are  found,  the  circulation  F,  of  the  flow  can  be  ob¬ 
tained  such  that 
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(52b) 


r. 


JL 

IT 


z 


jv(9,s)  SIM  H  &  J® 


with  V(0,s)- 


and 


For  an  example,  we  shall  take  a  slender  body,  the  sectional  form  of  which  is 
a  thin  symmetrical  Joukowski  aerofoil  shown  in  Fig.  1,  The  chord  length  and  the 
thickness  of  the  Joukowski  profile  are  represented  by  the  functions  c(»)  andcKs) 
of  5  ,  respectively.  We  shall  suppose  that  dhufcts)  is  a  constant,  and  that  the 

symmetrical  plane  of  the  body  has  a  small  angle  of  attack, ^  ,  against  the  uni¬ 
form  flow  (c.f.  Fig.  2).  Applying  the  general  expression  (52b)  to  this  problem, 
we  obtain 


(53) 


Since  both  the  edges  of  the  body  are  supposed  to  be  pointed, c(s)  vanishes  at  S-o 
and  S*|.  Remembering  Equation  (51e),  we  find  that  the  assvnnption  (50)  is  self- 
consistent  in  this  example. 


VIII  EXAMPLES 

Hereafter  we  shall  restrict  our  discussion  to  the  case  of  small  value  of  Rk. 
and  .  Expanding  the  right  hand  side  of  Equation  (27a)  into  the  R,i„  and  A  power 
series,  and  neglecting  the  order  of  square  of  Rn.  and  A  ,  we  have 

^=.B'(K*aC).o(RZ.&') 


Considering  the  character  of  the  Laplace  transform,  we  shall  obtain  the  linear  op¬ 
erators*/  and  luj  iva «  O  ^  in  Equation  (29a)  such  that. 


i.o(QT  R.4) 


(54a) 


(tty  '►C  )a,U)  =  a  i 


(54b) 
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Considering  the  character  of  the  Fourier  transform,  we  shall  obtain  the  operators 
and  such  that, 

^  (54c) 

103(1-0 (54d) 


In  the  subsonic  case,  both  the  operators 
transformation.  We  shall  obtain 


and 


S 


are  connected  with  the  Fourier 


(55a) 


C/03  )aJO 


W 


(55b) 


^  S  p  (55c) 

(55d) 

-B'o  i-H«' 

From  the  character  of  the  Fourier  transformation,  ^*(i)should  be  finite  for 
-  ^  s  S  .  Considering  Equation  (42b),  this  requirement  implies  that  2  <>) 

should  be  continuous  for-«»is§  ■♦■•o  .  Then  iV*)  should  vanish  at$-c  andS=J. 
Therefore  the  body  should  have  a  pointed  nose,  and  a  pointed  or  cylindrical  tail. 

We  shall  consider  two  examples. 

(i)  Slender  body  of  an  elliptic  cross  section  at  zero  incidence  (c.f.  Fig. 3). 
We  shall  suppose  that  no  singularity  is  at  the  axis.  The  sectional  form  of  the 
body  is  given  by 

_5l  .-i-  =  I 

(K4e)'  (n-er 

where  kli)  and  £(i)  are  the  functions  of  S  ,  and  f/K  is  supposed  to  be  a  small  con¬ 
stant.  Neglecting  the  terms  of  the  order  of  we  can  write  the  boundary 

conditions  in  the  following  form. 
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From  Equations  (42a)  and  (42b)  we  have 


«i(s)+tj(s) ^i(s)  =  -  t'-^k'e) 

From  Equation  (43a) ,  we  have 

-2(c.+  ^4o)-| 

=  -  z(^*  >►«*)»<.'*  K£. 

From  Equations  (43b)  and  (43e)  we  have 

/."a) 

K 

From  Equation  (43c)  we  have 


(56a) 

(56b) 

(56c) 

(56d) 

(56e) 

(56f) 


(i)  Vm 


%  (s)  4  ^  K*  -  t  (s)  K  €. 


(56g) 


From  Equation  (44)  we  have 


(56h) 


In  this  case  all  the  complex  functions  /t,  (s)  ,  etc.  are  reduced  to  the 

pure  real  functions  ,  i»(j)  etc.,  respectively.  Combining  these  eight  differ¬ 

ential  equations  (56a)-(56h)  with  the  seven  differential  equations  (41b),  (41c), 
(41e)  ,  (41f ) ,  (41g) ,  (41h)  and  (41 j),  we  can  determine  the  fifteen  arbitrary  func¬ 
tions  a,  ,  A,  ,  c,  ,  d,  ,  J,  .  ’  3«.  ’  ’  •^x  ,  .  and 
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Thus  we  obtain  the  solutions  ,  •-Ti-  ,  and  •■Qi*.  Remembering 

(31c)  and  neglecting  the  higher  order  terms,  we  have 

-t  U), 

where  4^,  *(  loj  4  2  'Z  (0  supersonic  case 


case 


where 


--0.3-t\(t*c)r2Vs)  for  the  subsonic 

=  JLl  Zi  (s)  i  S,  (s)  I 


K-Vdl. 

TV'S! 


for  the  supersonic  case 

'f,  >-cV<*hU  |3  are  replaced  by  'V  and  ^  for  the  subsonic  case 


with 


where 


(^1.:  I? 

|=/Z'(i)lo5k(5)ii 

-  f  (^0 

iK  ^  iii  /rzo 


Equation 

(57a) 

-X.J 

(57b) 


(57c) 

(57d) 

(57e) 


19 


The  velocity  field  in  the  sectional  plane  consists  of  the  diverging  or  converging 
radial  flow  and  the  quadrupole  like  flow.  The  perturbation  magnetic  field  is  pro¬ 
portional  to  R.^  .  The  outer  electric  field  is  represented  by  the  field  due  to  a 
quadrupole  at  the  origin.  The  lines  of  electric  force  for  the  inner  field  is  sim¬ 
ilar  to  the  stream  lines  of  the  fluid  flowing  along  the  inside  corner  of  right 
angle.  By  the  use  of  Equation  (45a),  the  surface  density  of  induced  charge  is 
represented  by 


2 


(58) 


The  surface  charge  takes  the  maximum  value  for  an  insulated  body  (R^"®).  As 
the  conductivity®^*  of  the  body  increases,  the  current  flowing  through  the  body 
increases.  As  a  result  the  charge  Ot  will  decrease.  Moreover  this  may  cause  the 
decrease  of  the  electric  field  given  by  Equations  (57d)  and  (57e).  Remembering 
Equations  (7a)  and  (7b),  we  can  introduce  the  electric  current  function  ♦  defined 
by 


1- 


(59) 


where 


The  outside  current  field  is  represented  by  superposing  a  quadrupole  field  on 
a  circular  field.  The  electric  field  is  shown  in  Fig.  4  schematically. 

(ii)  Body  of  revolution  at  a  small  incidence.  We  shall  consider  the  flow 
past  a  body  of  revolution  at  a  small  incidence, /tt  .  As  shown  in  Fig.  5,  the  x-axis 
is  taken  to  be  vertical.  By  the  use  of  the  similar  treatment  to  the  former  exampfe, 
we  can  solve  the  present  case.  Neglecting  the  order  of  square  of ,  4*  and  4^ 
are  approximately  represented  by  ^nd  ,  respectively.  The  functions 

and  are  obtained  by  replacing  Kli)  in  Equations  (57a)  and  (57b)  by  R/s)  . 
where  RCs)  denotes  the  radius  of  the  circular  cross  section.  The  functions 
and  y,  are  obtained  in  the  following  form 


(60a) 


with 


^ (ifv -  Ej!!. .  (Sf 


1  jl" 

U-o  *  J 

(60b) 


The  functions  ,X1.  ,si  ^ and  the  surface  charge  distribution  are  obtained 
in  the  following  form, 
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‘fc  = 

(60c) 

-- 

2  '  cT* 

Jl  = 

(60d) 

(60e) 

2  ^ 

=  - - /U  Sin  & 

(61) 

The  velocity  field  in  the  sectional  plane  consists  of  the  diverging  or  converging 
radial  field  and  the  dipole  like  field.  The  outside  electric  field  is  that  due  to 
a  dipole  at  the  origin,  and  the  inside  electric  field  is  uniform.  The  coefficient 
has  the  same  physical  meaning  as  in  the  previous  example.  The  electric 
field  and  the  surface  charge  distribution  are  schematically  shown  in  Fig.  5. 


IX  TOTAL  FORCE 


The  force  exerted  on  a  body  in  the  uniform  flow,  and  the  applied  magnetic 
field  consists  of  the  hydrodynamic  pressure  and  the  electro-magnetic  force.  Since 
there  is  no  applied  electrostatic  field  in  the  present  problem,  the  total  electro¬ 
static  force  exerted  on  the  body  should  vanish.  As  the  Lorentz  body  force  exerts 
on  the  induced  current  in  the  interior  of  the  body  by  the  applied  magnetic  field, 
it  should  vanish  for  an  insulated  body.  According  to  the  electromagnetics,  the 
Lorentz  body  force  can  be  expressed  by  the  surface  integral  of  the  Maxwell  stress 
tensor.  By  the  use  of  the  momentum  theorem,  the  total  force JF  can  be  expressed  by 
the  integral  over  a  control  surface,  5,  say,  such  that 


where  y/  denotes  the  velocity  vector 
vector  normal  to  the  control  surface 


,  (62) 

s. 


denotes 


the  outward  directing  unit 


(A)  DRAG 

The  drag  0  is  the  s-component  of  the  total  force  .  After  the  similar 
treatment  to  the  conventional  case  ^  s-component  of  Equation  (62)  can  be 
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reduced  to 


-,2CLj'[4(H»Il)(||)jr]ds 

•  Cli) 

rV  rdcp  (\3<f  I3V  "S^C^  /3(P  . 

-zaj  4L»^J 

•  ccO 

-  ^  [l7,'(dVds»ofa‘R-) 


5.r  M.  >  1 


with 


if.*  U.“  LV 


(63) 


where  C(J)  denotes  the  cross  sectional  form,  and  H  and  denote  the  outward 
directing  normal  and  the  counter-clockwise  vectors  of  C,  respectively.  Since  the 
drag  given  by  Equation  (63)  corresponds  to  the  integral  of  the  hydrodynamic  pres¬ 
sure  over  the  surface,  the  Lorentz  force  does  not  make  any  contribution  to  the 
total  drag  within  the  accuracy  of  the  present  approximation.  The  first  term  of 
the  right  hand  side  of  Equation  (63)  represents  the  conventional  supersonic  drag. 
In  the  subsonic  case,  this  term  is  dropped,  and  other  terms  do  not  change.  For 
example  we  shall  calculate  the  drag  for  a  slender  body  of  revolution  at  a  small 
incidence.  Substituting  the  solutions  (60a),  (60b)  and  (60d)  into  Equation  (63), 
we  have 


_  a 


TT 


I 

5-5 


'  z 


•A 


Ho  >  I 


(64a) 
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The  first  term  of  the  right  hand  side  of  Equation  (64a)  represents  the  conventional 
supersonic  drag,  which  is  the  order  of  0(2',)^  .  The  sign  of  the  second  term  cor¬ 

responds  to  the  s-distribution  of  the  cross  sectional  area,  and  the  order  of  this 
term  is  OiAX*)  .  The  third  term  represents  the  mutual  effect  between  the  com¬ 
pressibility  and  the  conductivity.  This  term  always  takes  the  positive  sign  and 
is  the  order  of  The  last  term  represents  the  so-called  induced  drag, 

which  is  the  order  of  o('ai'2r*)  .  The  induced  drag  takes  the  negative  sign  in  the 
present  case.  Since  the  third  term  predominates  among  the  last  three  terms  of 
Equation  (64a),  the  total  drag  increases  by  the  electromagnetic  effect.  Putting 
in  Equation  (64a),  we  shall  obtain  the  drag  for  a  body  of  revolution. 

This  expression  coincides  with  the  result  due  to  Ando^^2)^  except  only  one  term. 

The  difference  results  from  his  treatment,  where  he  omitted  the  last  integral  in 
the  pressure  formula  (8a)  before  substituting  it  into  Equation  (62).  In  the  con¬ 
ventional  case,  the  term  due  to  the  compressible  effect  vanishes  unless  the  area 
of  the  base  section  2  0)  does  not  vanish.  In  the  present  case,  the  terms  due  to 
the  electromagnetic  effect  still  remains  finite  for  a  body  with  a  pointed  tail. 


In  the  subsonic  case,  the  drag  is  represented  by 


(64b) 


where  3  j- 


F«r  Mo  I 


The  conventional  drag  vanishes.  The  behavior  of  other  terms  are  similar  to  those 
of  the  supersonic  case. 

(B)  LATERAL  FORCE 


^  By  the  use  of  Equation  (62),  the  x-and  the  y-components  of  lateral  force, 
^  ,  Y*  ,  are  calculated  such  that 

X  -t  -i  Y  « 


with 


+  0 


(  a',  K.*) 


(65a) 


X  = 


ur  LV 


where  ^  denotes  the  angle  between  the  normal  vector  “Vj  of  the  curve  C  and  the 
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x-axis. 


Applying  Equation  (65a)  to  the  body  of  revolution  at  a  small  incidence, 
we  obtain 


Y-  o 


(65b) 


The  equation  (65b)  implies  that  the  lift  takes  the  negative  sign  for  a  positive 
value  of  the  incidence,  /ul.  .  This  negative  lift  corresponds  to  the  negative  in¬ 
duced  drag  discussed  above.  For  the  sake  of  considering  the  physical  meaning  of 
the  negative  lift,  we  shall  decompose  the  lift  to  the  two  integrals  of  the  pres¬ 
sure  and  the  Maxwell  stress  tensor  over  the  surface  of  the  body.  The  integral  of 
the  pressure  is  f.ven  by 


The  integral  of  Maxwell  stress  tensor  is  given  by 

-2  S  f  [ 

First  we  shall  consider  an  insulated  body.  We  shall  suppose  is  positive  for 
O  <  S  c  9,  ,  and  negative  for  S,  <  S  < «  .  For  the  region,  where  'ZVr)is  positive, 

the  velocity  field  in  the  sectional  plane  consists  of  the  dipole  like  and  the  di¬ 

verging  radial  fields  (c.f.  Fig,  6).  As  a  result,  the  induced  current  flows 
clockwise  around  the  circular  section  C,  The  Lorentz  force  due  to  the  correlatbn 
between  the  induced  current  and  the  uniform  applied  magnetic  field  acts  in  the 
inward  direction  normal  to  the  circle.  Remembering  the  velocity  field  shown  in 
Fig,  6,  it  is  realized  that  the  velocity  due  to  the  dipole  and  the  Lorentz  force 
take  the  similar  directions  with  each  other  on  the  upside  of  the  body,  but  they 
take  the  opposite  directions  on  the  downside.  We  shall  consider  a  couple  of  small 
volume  elements  of  fluid  flowing  along  stream  lines  of  the  upside  and  the  down¬ 
side  of  the  body.  The  work  performed  by  the  Lorentz  force  when  the  fluid  passes 
along  the  stream  line  of  the  upside  from  o  to  3  ,  should  be  larger  than 

the  work  performed  when  the  fluid  passes  along  the  stream  line  of  the  downside. 

Thus  the  pressure  on  the  upside  is  larger  than  that  on  the  downside  for  the  same 
section,  S  .  After  the  position  S,  ,  2  Vv  is  negative.  The  velocity  field  in 
the  sectional  plane  consists  of  the  dipole  like  field  and  the  converging  radial 
field.  Then  the  induced  current  flows  counterclockwise,  and  the  Lorentz  force 
acts  in  the  outward  direction.  As  a  result,  the  work  performed  on  the  upside  is 
smaller  than  the  work  performed  on  the  downside.  Thus  the  excess  pressure  on  the 
upside  decreases  as  S  increases  ( 5,  <  i  ) ,  and  the  pressures  of  the  both  sides 
take  the  same  value  at  the  tail,  3=(  .  Since  the  upside  pressure  is  larger  than 
the  downside  pressure,  the  integral  of  the  pressure  distribution  takes  the  nega¬ 
tive  sign  (c.f.  Equation  (66a)). 


(66a) 


(66b) 


Next  we  shall  consider  a  conducting  body.  Remembering  Equation  (60e),  we  have  the 
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inside  current  density  such  as 


The  Lorentz  force  due  to  the  correlation  between  the  uniform  horizontal  in¬ 
duced  current  and  the  applied  magnetic  field  acts  downward  in  the  vertical  direction 
(c.f.  Equation  (66b)).  Thus,  both  the  integrals  of  the  pressure  and  the  Maxwell 
stress  tensor  contribute  to  the  negative  lift. 
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Figure  2 


Meridian  section  of  slender  body  at  small  incidence 
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Figure  5  Induced  electric  field  and  surface  charge 

distribution  for  body  of  revolution  at  small 
incidence 
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